REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  foj»r  rawing  instructions,  scan 
completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  inform* 

Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlin 

other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valio  - - 

ABOVE  ADDRESS.  _ 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

04/30/01  Final  Report  _ 


4.  TITLE  AND  SUBTITLE 

Ultrashort  Pulse,  Monolithic  Modelocked  Lasers  for  WDM  Systems 


AFRL-SR-BL-TR-01- 

05u,-2^ 


3.  DATES  COVERED  (From  -  To) 
07/15/96  -  03/14/00 


5a.  CONTRACT  NUMBER 
F49620-96-1-0371 
5b.  GRANT  NUMBER 


6.  AUTHOR(S) 

P.  Daniel  Dapkus 


5C.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Electrical  Engineering  -  Electrophvsics 
University  of  Southern  California 
Los  Angeles,  CA  90089-0243 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Research 
801  North  Randolph  Street,  Room  732 
Arlington,  VA  22203-1977 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release:  Distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

A  program  of  research  to  investigate  the  use  of  modelocked  lasers  as  sources  for  WDM  systems  is 
reported.  The  objective  of  the  program  was  to  demonstrate  modelocked  lasers  with  ultrashort  pulse 
widths  that  contain  as  many  as  16  phase  locked  longitudinal  modes.  Analysis  was  carried  to  understand 
the  role  of  self  phase  modulation  on  the  spectral  range  of  modes  that  could  be  phase  locked  in  such  a 
device.  The  results  of  the  analysis  suggest  that  this  fundamental  limitation  must  be  addressed  by 
engineering  the  round  trip  phase  dispersion  over  a  broad  spectral  range  to  achieve  a  broad  enough 
spectral  width  to  successfully  fabricate  monolithic  WDM  sources  that  will  compete  with  current 
technology.  Experimental  monolithic  passively  modelocked  colliding  pulse  modelocked  lasers  were 
constructed  during  the  course  of  the  program.  Devices  with  1.2  psec  pulse  width  operating  with  a  pulse 
repetition  rate  of  -  40  GHZ  were  demonstrated.  These  devices  had  longitudinal  mode  spectra  containing 
up  to  16  phase  locked  modes  but  the  amplitude  variation  over  the  mode  spectrum  precluded  their  use  as 
independent  WDM  spectral  sources. 


15.  SUBJECT  TERMS 

Optical  Communications,  Photonics,  Wavelength  Division  Multiplexing,  Modelocked  Lasers,  Semiconductor 
Lasers 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

P .  Daniel  Dapkus 

a  REp0RT 

Unclassif ied 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

UL 

28 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

(213)  740-4414 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


20010623  122 


Introduction 

This  report  describes  the  result  of  an  experimental  and  theoretical  study  of  the  use  of  mode  locked  lasers  as 
multiple  wavelength  sources  for  WDM  systems.  The  effort  involved  modeling  of  mode  locked 
semiconductor  lasers  with  special  emphasis  on  understanding  the  role  of  self  phase  modulation  on  the 
locking  stability  and  range  of  passively  mode  locked  lasers.  Our  objective  in  undertaking  this  work  was  to 
extend  the  spectral  range  of  mode  locked  lasers  so  that  several  modes  ( — 16)  could  be  locked  to  produce  a 
manifold  of  equally  spaced  laser  wavelengths  to  be  used  in  a  WDM  system. 

Progress  in  this  program  was  severely  hampered  by  the  untimely  illness  of  the  graduate  student  involved  in 
the  program.  Approximately  1.5-2  years  into  the  program,  the  student  who  was  carrying  out  the 
experimental  device  work  suffered  a  debilitating  illness  that  ended  his  graduate  work  and  severely  truncated 
the  experimental  effort.  This  prevented  us  from  achieving  some  of  the  program  objectives  listed  below. 
However,  we  were  able  to  accomplish  some  of  the  program  objectives  and  to  define  to  likelihood  of  success 
for  the  approach  originally  proposed. 

Mode  locking  is  effected  in  a  laser  by  modulating  the  net  gain  in  a  laser  at  a  frequency  equal  to  the 
frequency  separation  of  the  longitudinal  modes  of  the  laser.  This  has  the  effect  coupling  nearest  neighbor  or 
next  nearest  neighbor  modes  together  and  locking  the  phase  of  these  modes  to  one  another.  The  net  effect  of 
this  phase  locked  system  of  modes  is  that  they  behave  as  a  collective  ensemble  of  photons  that  reaches 
threshold  at  a  value  lower  than  would  occur  for  a  conventional  single  mode  operation.  The  coherent  phase 
locked  ensemble  of  modes  in  the  frequency  domain  leads  to  a  repetitive  short  pulse  operation  in  the  time 
domain.  The  repetition  frequency  of  the  pulses  is  equal  to  the  round  trip  frequency  of  the  photons  in  the 
cavity,  c/(2nL).  In  principle  such  a  laser  can  be  used  as  a  high  repetition  rate  short  pulse  generator  if  all  of 
the  coherent  modes  are  used  in  the  output  pulse  or  if  the  individual  modes  are  filtered  separately  from  the 
output  as  a  CW  multi-wavelength  source.  In  either  form,  this  laser  is  useful  as  a  signal  generation  source  for 
photonic  systems.  In  principle,  because  the  gain  spectrum  of  a  semiconductor  laser  is  so  wide,  the  device 
should  be  able  to  exhibit  many  modes  (>  1 6)  and  ultrashort  pulses  (<0.1  ps).  In  practice  neither  of  these 
limits  has  been  approached.  The  key  reason  was  believed  to  be  dispersion  in  the  index  of  refraction  in  a 
semiconductor  laser.  Such  dispersion  leads  to  unequal  mode  separation  and  ineffective  locking. 

Prior  to  the  beginning  of  this  program  there  had  been  studies  of  the  mode  locking  behavior  of  semiconductor 
lasers  to  exploit  their  beneficial  properties.  The  best  results  achieved  had  been  those  reported  in  [1]  with  a 
colliding  pulse  mode  locked  (CPM)  laser.  In  these  devices,  two  counter  propagating  pulses  collide  in  a 
center  saturable  absorber  to  saturate  the  loss  in  this  device.  The  presence  of  the  two  pulses  increases  the 
power  available  to  saturate  the  absorber  on  the  leading  edge  of  the  pulses  and  thus  to  effectively  shorten  the 
pulses.  In  this  work  a  pulse  width  of  ~  1  ps  was  achieved  with  approximately  8  locked  modes. 

The  objective  of  this  program  was  to  improve  on  these  results  by  designing  a  cavity  with  minimal  round  trip 
dispersion  and  with  minimal  dynamic  phase  dispersion.  This  was  to  be  accomplished  by  varying  the  gain 
peak  and  thus  the  index  dispersion  on  a  local  basis  along  the  cavity  by  using  a  spatially  varying  energy  gap 
in  the  active  region.  It  is  believed  that  such  a  gain  region  can  be  grown  using  selective  area  epitaxy  and  that 
the  cavity  can  be  designed  to  minimize  static  and  dynamic  phase  dispersion.  To  accomplish  this  objective  a 
program  in  which  extensive  theoretical  modeling  was  to  be  undertaken  to  guide  the  device  design  and  the 
experimental  research.  Furthermore  a  program  of  experimental  research  in  which  studies  of  selective 
epitaxial  growth  to  form  the  spatially  varying  active  region  were  undertaken  to  augment  device  design, 
fabrication  and  evaluation  experiments. 
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2.0  Experimental  Studies 


2.1  Laser  design 

The  basic  laser  design  to  be  used  in  this  study  is  the  colliding  pulse  mode  locked  (CPM)  laser 
shown  schematically  in  Figure  1.  It  consists  of  a  gain  region,  G,  a  modulator  region,  M,  and  a 
saturable  absorber  region,  SA.  These  regions  may  be  sequentially  placed  along  the  cavity  or 
positioned  symmetrically  about  the  centered  SA  as  shown  in  the  diagram. 


Figure  1  Schematic  diagram  of  CPM  laser  to  studied  in  this  program.  G  are  the  gain  sections,  SA  the  saturable 
absorber,  and  M  the  modulators. 


Symmetric  placement  is  preferred  to  produce  the  shortest  pulse  performance.  In  this  device  design 
the  gain  regions  provide  amplification  for  the  optical  wave,  the  modulator  sections  provide  a  time 
dependent  loss  the  frequency  of  which  is  matched  to  the  round  trip  frequency  of  the  light  wave, 
c/2nL,  and  the  SA  region  provides  an  intensity  dependent  loss  mechanism  to  sharpen  the  leading 
redge  of  a  propagating  pulse.  The  waveguide  of  the  device  is  formed  either  by  forming  a  buried 
heterostructure  (as  schematically  illustrated)  or  through  the  use  of  a  ridge  waveguide.  Individual 
contact  is  made  to  each  of  the  regions  of  the  device  and  the  contacts  common  are  inteconnected  to 
the  same  contact  pad.  The  contact  pads  are  isolated  from  the  substrate  by  an  insulating  SiNx  layer. 

The  specific  lengths  of  the  individual  regions  was  to  be  determined  by  experiments  and  analytical 
modeling.  The  total  length  of  the  device  was  chosen  to  coincide  with  a  a  round  trip  frquency  of  40 
GHz.  This  resulted  in  a  longitudinal  mode  separation  of  80  GHz  since  a  symmetric  CPM  device 
operates  with  a  mode  separation  twice  that  is  twice  the  round  trip  frequency.  40  GHz  was  also  the 
limit  of  available  RF  generators  for  the  modulator  drive. 

The  active  region  of  the  device  is  formed  by  growth  of  strained  multiple  quantum  well  (QW)  layers 
designed  to  emit  at  1.55  pm.  A  significant  part  of  the  program  was  designed  to  exploit  selective 
area  growth  (SAG)  to  be  able  to  tailor  the  active  regions  to  optimize  device  performance.  This  will 
be  described  in  detail  in  the  next  section.  However,  the  basic  principle  involves  varying  the  QW 
thickness  on  a  spatial  scale  sufficient  to  vary  the  effective  energy  gap  of  the  device  regions  to 
optimize  their  function  in  the  device  operation.  This  is  accomplished  by  patterning  a  dielectric 
layer  placed  on  the  wafer  prior  to  growth  of  the  active  region.  The  pattern  is  chosen  to  mask  the 
region  around  the  areas  in  which  growth  is  to  occur.  Owing  to  suppression  of  the  growth  in  the 
mask  regions  some  of  the  materials  that  would  have  deposited  there  transport  in  the  vapor  and  on 
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the  surface  to  enhance  the  growth  rate  in  the  openings.  The  growth  rate  enhancement  is  determined 
by  the  relative  dimensions  of  the  open  and  masked  areas  and  leads  to  thicker  quantum  wells  with 
smaller  effective  energy  gaps.  By  tailoring  the  pattern,  we  can  then  control  the  effecive  bandgap  in 
the  active  region  spatially.  For  example,  It  may  be  desirable  to  increase  the  energy  gap  of  the 
saturable  absorber  to  allow  it  to  be  operated  at  higher  voltage  for  rapid  sweep  out  of  the  carriers. 
Similar  considerations  may  apply  to  the  modualtor.  The  next  section  describes  the  work 
accomplished  in  this  program  to  allow  us  to  selectively  grow  the  various  regions  of  the  device. 

The  fabrication  of  the  device  used  standard  procedures 
except  for  controlling  the  length  of  the  device.  Ridge 
waveguides  were  fabricated  by  dry  etching  a  ridge  to 
define  the  waveguide  down  to  an  etch  stop  layer  above 
the  guide  layers  of  the  structure.  The  structure  was 
covered  by  SiNx  and  openings  for  the  contacts  to  the 
stripes  were  defined  and  etched.  Ti/Pt?Au  ohmic 
contacts  were  applied  to  the  surface  through  a 
photoresist  mask  and  the  residual  metal  was  lifted  off  to 
define  all  of  the  metal  contacts  to  the  various  reions 
simultaneously.  Openings  to  define  position  of 
subsequent  cleaves  were  then  mad  in  the  wafer.  These 
openings  were  narrow  stripes  aligned  perpendicualr  to 
the  cavity  at  the  prsecribed  separation  to  yield  a  1  mm  device. 
An  etchant  that  defines  (111)  planes  during  etching  was  used 
to  form  sharp  V-grooves  at  the  position  of  the  openings.  Thes 
V-grooves  act  as  stress  risers  during  the  subsequent  cleaving 
operation  and  yield  devices  lithographically  controlled 
lengths.  The  wafer  was  then  thinned  from  the  back  to  ~  100 
pm  and  n-type  ohmic  contacts  were  applied  to  the  back  of  the 
device.  At  this  time  the  wafer  was  cleaved  bars  by  applying 
bending  stress  to  the  wafer.  Each  bar  contains  several  devices 
that  can  be  subsequent  sawn  from  the  bar  or  probe-tested. 
Figure  2  shows  a  micrograph  of  three  lasers  prior  to  cleaving. 
These  lasers  are  designed  to  be  passively  modelocked  devices  and  thus  have  no  modulator  sections. 
The  saturable  absorber  is  40  pm  long.  Figure  3  shows  the  end  of  a  bar  of  devices  that  have  been 
cleaved  using  the  V-groove  location  control  technique. 


Figure  2  Top  view  od  CPM  Modelocked 
lasers  made  in  this  program 


Figure  3  View  of  edge  of  laser  bar 
showing  precision  cleave  grooves 


2.2  Materials  Growth 

2.2.1  Calibration  of  growth  with  TBAs  and  TBP 

Significant  progress  towards  the  required  selective  area  materials  technology  to  fabricate  lasers 
with  spatially  varying  energy  gap  along  the  active  region  was  made  during  the  course  of  this 
program.  Selective  area  growth  involves  the  patterning  of  the  wafer  with  an  amorphous  dielectric 
to  inhibit  growth  in  certain  areas.  The  growth  in  adjacent  areas  is  then  enhanced  because  they 
receive  a  greater  supply  of  reactants.  Inhibition  of  growth  in  the  masked  areas  cause  reactants  that 
would  otherwise  deposit  on  the  masked  areas  to  diffuse  laterally  in  the  gas  phase  or  migrate  on  the 
surface  to  be  incorporated  in  the  unmasked  areas.  This  is  shown  schematically  in  Fig.4.  The  mask 
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dimension,  wm,  relative  to  the  open  area,  w0, 


(1)  Vertical  vapor  diffusion 

(2)  Lateral  vapor  diffusion 

(3)  Surface  migration 


Figure  4  Schematic  diagram  showing  the  mechanisms  operating  in  SAG  to  produce  enhanced 
growth  rate  on  the  mask  opening. 
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Figure  5  Typical  bubbler  design  for  MOCVD  sources  with  multiple 
outputs. 
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Figure  6  TBP  bubbler  pickup  efficieny  vs.  flow  rates 


determines  the  enhancement  of  the 
growth  rate  in  the  open  area  over 
that  which  would  be  observed  for  an 
unmasked  substrate.  Because  the 
various  constituents  of  the  lattice 
diffuse  and  migrate  at  different 
rates,  there  is  also  the  possibility  of 
a  composition  change  that 
accompanies  this  selective  growth. 
To  understand  the  degree  to  which 
such  effects  occur  it  is  important  to 
carefully  calibrate  the  flux  of 
reactants  in  the  system.  This  is 
particularly  important  for  our 
reactor  since  we  are  using 
organometallic  reactants  for  all  of  the 
sources,  trimethylindium  (TMIn), 
trimethylgallium  (TMGa), 
tertiarybutylarsine  (TBAs)  and 
teriarybutylphosphine  (TBP).  A 
careful  study  of  the  transport  of  TBP 
in  arsine  was  carried  out  because  the 
growth  system  involved  a  dilution 
system  to  be  able  control  the  TBP 
flux  over  a  wide  range.  This  portion 
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ITBP1  =  5.63  seem 


d/C  «  =  4.43 


of  the  system  is  shown  schematically  in  Fig  5. 
l  _  l  c^p  Passing  a  carrier  gas,  [In],  through  the  TBP  liquid,  the 

crap  “  [TBP]  [  c'h’  saturated  carier  is  diluted  with  another  flow  of  carrier 

7|  gas,  [Dil],  to  reduce  and  control  the  concentration  of 
15  "  Itmci  =  o.o2o  seem  the  TBP  in  the  carrier  gas.  This  diluted  gas  stream  is 

o  °82 v4=oC5o  s'  then  sampled  by  the  Epison  to  measure  the 

.2  io  -  yt  concentration  of  TBP  in  the  gas  and  the  output  is  split 

2  /  between  two  out  put  mass  flow  controllers  ,  [Outl] 

g  5  -  /itbpi  =  5.63  seem  and  [Out  2].  Since  these  mass  flow  controllers  were 

g  not  specifically  calibrated  for  TBP  and  detailed 

0  -  CiBp/c  H  _  4.43  calibration  is  necessary  to  determine  the  precise  flow 

/C  /  ’  of  TBP  that  enters  the  reactor.  The  pickup  from  the 

'5  ^ ^ ^ ^ ^  ^  bubblers  was  measured  by  determining  the 

Apparent  MFC  setting  (seem)  concentration  of  reactant  in  the  gas  flow  from  the 

bubbler  for  flows  and  temperatures.  This  data  shown 
in  Fig.  6  indicates  that  pickup  of  TBP  is  linear  with 
e'fson  7  Calibratlon  of  TBP  bubbler  witb  the  dilution  flow  over  the  useful  range.  The  output 
epison‘  mass  flow  controllers  were  calibrated  by  growing  a 

series  of  InGaAsP  samples  of  the  same  composition  under  different  TBP  flow  conditions.  Since  all 
other  variables  are  fixed,  the  amount  of  TBP  entering  the  reactor  was  known  to  be  the  same  for  all 
samples.  The  data  and  appropriate  calibration  equations  are  shown  in  Fig.  7.  This  data  can  only  be 
interpreted  if  the  specific  heat  of  the  TBP  is  15.3  R.  Theoretical  calculations  also  predict  the  heat 
capacity  to  be  15.5R,  an  excellent  agreement.  These  experiments  are  crucial  to  our  knowing  the 
t.o _  absolute  mole  fraction  input  of  TBP  and  TBAs. 


Apparent  MFC  setting  (seem) 

Figure  7  Calibration  of  TBP  bubbler  with 
epison. 
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The  absolute  mole  fractions  of  the  reactants  are  very 
important  variables  to  develop  an  understanding  of  the 
selective  area  growth  enhancement  on  patterned 
substrates  and  to  predict  the  conditions  to  grow 
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Figure  8  Calibration  of  the  In 
composition  of  the  Q  material  vs.  gas 
phase  composition. 

appropriate  structures  on  planar  substrates  as 
well.  For  example,  we  expect  the  column  III 
composition  to  be  linearly  dependent  on  the 
gas  phase  reactant  ratio.  This  trend  is  clearly 
shown  in  Figure  8,  where  the  In  composition 
and  the  As  composition  of  quaternary 
materials  are  plotted  as  a  function  of  the  gas 


As/P-Efficiency  = 


(As/P)Sl 

(As/P)v. 


[TBP 

[TBA 


ITBPl/dTMN  +  ITMCI) 

Figure  9  As/P  efficiency  in  the  solid  vs.  V/III 
ratio. 
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phase  composition  and  the  growth  temperature  at  a  fixed  gas  phase  composition,  respectively.  The 
close  correlation  of  the  data  with  the  expected  behavior  results  because  of  the  accurate  calibration 
we  achieved  in  these  studies. 

The  column  V  composition,  however,  is  a  much  more  complicated  function  of  the  growth 
parameters.  In  Figure  9,  the  V  incorporation  efficiency  is  shown  as  a  function  of  the  TBP  /  III  gas 
phase  ratio.  The  data  show  that  the  As  incorporates  more  efficiently  at  higher  TBP  mole  fractions 
independent  of  the  growth  rate  for  either  1 .2  pm  or  1 .5  pm  materials.  This  makes  control  of  the 
composition,  and  thus  strain,  a  complicated  function  of  the  growth  parameters. 


2.2.2  Quantum  Well  Growth 

To  achieve  reproducible  growth  of  quantum  well  active  regions,  we  have  chosen  a  design  that 
requires  the  use  of  two  Q  compositions  and  InP.  This  design  is  shown  in  Figure  10.  The  quantum 
wells  consist  of  Q(A=1.7  pm)  wells  and  Q  (X,  =  1.2  pm)  barriers.  The  1.5  pm  emission  is  achieved 
by  controlling  the  quantum  size  effect  through  the  well  thickness.  The  Q  (A,  =  1.2  pm)  barriers 
insure  that  a  high  confinement  potential  exists  to  confine  the  carriers  at  room  temperature. 


InP - 

0(  1.2pm) 

well 


Figure  10  Schematic  diagram  of  MQW  structure  used  in 
this  program. 


High  quality  QW  structures  have  been  grown  with  this  design  and  characterized  to  calibrate  the 
emission  wavelength  vs.  well  thickness.  The  photoluminescence  of  three  different  samples  with 
different  well  widths  are  shown  in  Figure  11.  The  spectra  shown  are  characteristic  of  high  quality 
materials.  The  spectra  also  show  the  variation  of  the  emission  wavelength  as  a  function  of  the  well 
width.  This  is  summarized  in  part  b.)  of  the  figure. 


Wavelength  (pm) 


a.) 

Figure  1 1  a)  PL  spectra  of  QW  samples  with  different 
width  for  strained  QW  layers  with  1 .2  pm  barriers. 


width.  B)  Plot  of  emission  wavelength  vs  well 


The  quality  of  the  materials  is  further  verified  by  the  excellent  correlation  of  the  X-Ray  diffraction 
spectra  measured  on  these  structures  as  shown  in  Figure  12.  This  figure  shows  the  measured  and 
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simulated  x-ray  diffraction  spectra  from  the  MQW 
structures  whose  PL  is  shown  in  Figure  11.  The 
presence  of  high  order  diffraction  peaks  and  the 
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Figure  12  Xray  rocking  curve  of  Laser 
structure. 

excellent  correlations  of  peak  position  and 
relative  intensity  attest  to  the  quality  of  the 
materials.  Figure  13  shows  the  PL  spectrum  of  a  MQW  sample  grown  with  the  appropriate 
structure  to  emit  at  1.55  pm.  The  sample  exhibits  a  very  narrow  linewidth  indicative  of  high  quality 
material.  Based  on  these  characterizations  it  was  expected  that  the  materials  were  sufficiently  high 
quality  to  support  the  fabrication  of  high  performance  devices. 


Figure  13  PL  spectra  of  MQW  sample  showing 
narrow  linewidth 


2.2.3  Buried  Heterostructure  Device  Fabrication 

Since  buried  heterostructure  (BH)  lasers  are  known  to  exhibit  the  lowest  threshold  currents  among 
the  many  designs  available  for  implementation  as  a  mode  locked  laser  and  since  the  structures  are 
most  suitable  for  implementing  single  lateral  mode  designs,  effort  was  devoted  to  development  of 

these  structures  for  our  ML  lasers.  The  key  task 
in  fabricating  a  high  quality  BH  structure  is  to 
regrow  a  current  blocking  layer  over  the  mesa 
to  be  etched  in  the  active  region  that  defines  the 
waveguide.  An  example  of  the  mesa  ridge  to  be 
buried  is  shown  in  Figure  14.  The  mesa  is  deep 
enough  to  penetrate  the  active  region  of  the 
device  and  narrow  enough  to  ensure  that  the 
resultant  mode  will  be  single  mode  after  the 
mesa  is  buried  in  materials  with  comparable 
index  of  refraction.  In  addition,  the  burying 
layers  must  block  current  flow  so  that  all  of  the 
current  flows  through  the  active  region.  To 
accomplish  this  a  series  of  alternately  doped 
layers  are  grown  selectively  around  the  mesa.  During  the  growth  of  the  current  blocking  layer 
(CBL),  the  top  of  the  mesa  is  covered  with  SiNx  to  prevent  growth  on  the  top  oif  the  mesa.  This  can 
result  in  complicated  growth  morphology  and  current  leakage  paths  around  the  mesa.  The  desired 
growth  morphology  is  shown  in  Figure  15a.  Here  we  illustrate  the  importance  of  having  the  first 
grown  p  layer  discontinuous  with  the  p  layer  n  the  top  of  the  mesa  to  minimize  the  possibility  of  a 
current  leakage  path.  The  two  CBL  structures  shown  in  the  microphotographs  of  b.)  and  c.)  shown 


Figure  14  SEM  Image  of  as-etched  BH  mesa 
before  burying. 
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structures  with  the  N  layer  of  the  CBL  at  0.2  jam  and  0.6  (am  thickness,  respectively.  In  the  latter 


c. 


Figure  1 5  a)  Schematic  diagram  of  CBL  structure.  B) 
and  c)  SEM  cross  ectional  images  of  two  BH  laser 
structures 

device  with  the  thicker  n  layer  in  the  CBL.  These 
results  suggest  the  use  of  thinner  n  -  layers  in  the 


case  note  the  formation  of  large  growth  facets  in 
the  n-layer  around  the  mesa.  These  facets 
disrupt  the  subsequent  over  growth  of  the  top 
contact  layer  and  apparently  create  a  current 
path  that  influences  the  leakage  currents  of  the 
device  as  shown  in  Figure  16.  Here  we  show  the 
I  V  characteristics  of  these  two  BH  devices. 
Note  the  large  leakage  in  the  reverse  direction 
and  the  lower  forward  turn  on  voltage  of  the 


Figure  16  I-V  characteristics  of  BH  laser  structures 
showing  effect  of  thick  n  CBL  layer. 

CBL. 


Optimized  CBL  layers  and  device  structure  as  well  as  the  resultant  I-V  characteristics  are  shown  in 


Voltage( 

a.)  b.) 


Figure  17  a)  Cross  section  of  BH  laser  structure.  B)  I  -  V  characteristics  of  two  devices  showing  high  leakage 
in  device  with  CBL  grown  at  a  high  rate. 
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Figure  17a.)  and  b.).  The  data  indicate  that  as  long  as  the  n-layer  of  the  CBL  is  grown  at  low  rates 
so  that  its  thickness  near  the  edge  of  the  mesa  can  be  controlled,  low  leakage  structures  can  be 
grown.  These  devices  result  in  low  threshold  room  temperature  CW  operation.  The  L-I 
characteristics  and  spectra  of  devices  operating  under  CW  conditions  are  shown  in  Figure  18  a)  and 
b).  These  devices  have  2  pm  wide  active  regions  with  four  QW’s.  They  operated  CW  with  an 
emission  wavelength  near  1.55  pm. 


1.545  1.550  1.555  1.560 

wavelength  (//m) 


a)  b) 

Figure  18  a)  Light  output  versus  current  for  BH  lasers  of  different  length,  b)  Spectra  of  the  same  lasers  operating 
under  CW  conditions. 


2.2.4  Selective  Area  Growth 

As  discussed  in  the  Introduction  and  in  other  sections,  the  program  plan  called  for  us  to  explore  the 
benefits  of  a  spatially  varying  active  region  energy  gap  on  the  performance  of  ML  lasers.  As  a 
result  significant  experimentation  was  undertaken  in  the  program  to  develop  techniques  for  varying 
the  energy  gap  of  the  active  region  along  the  active  stripe  use  of  selective  area  growth  techniques. 
Although  they  were  not  employed  in  real  ML  lasers  these  techniques  are  central  to  many  integrated 
device  concepts  and  the  results  will  be  applied  in  many  areas. 

The  basic  concept  of  selected  area  growth  (SAG)  is  illustrated  in  Figure  4.  A  dielectric  mask 
inhibits  the  deposition  of  material  in  certain  areas.  This  increases  the  growth  rate  in  adjacent  areas 
because  the  gas  phase  reactants  diffuse  to  the  open  areas  in  the  mask  driven  by  a  concentration 
gradient  and  because  reactants  that  land  on  the  mask  migrate  over  the  mask  area  to  the  open 
regions.  Gas  phase  diffusion  leads  to  a  more  or  less  uniform  increase  in  the  deposition  rate  over  the 
entire  open  area,  but  the  surface  migration  leads  to  non  uniform  increases  at  the  margins  of  the 
open  areas.  Owing  to  surface  kinetic  phenomena,  slow  growth  planes  may  be  identified  in  the 
profile  of  the  deposit  in  the  opening.  In  Figure  19  we  show  cross  sectional  micrographs  of  an  InP 
deposits  on  an  opening  in  a  mask  on  an  InP  substrate.  Periodic  ultrathin  layers  of  InGaAs  are 
deposited  to  decorate  the  growth  habit.  The  two  cross  sections  are  markedly  different  because  the 
opening  edges  are  aligned  along  different  crystalline  directions.  In  19  a)  the  stripe  openings  are 
aligned  along  the  <1 1 0>  direction.  Under  these  circumstances,  slow  growth  (1 1 1  A)  planes  are 
formed  at  the  edges  that  define  the  sidewalls  of  the  mesa.  In  the  orthogonal  direction  (19  b)),  no 
such  plane  is  present  and  the  crystal  grows  laterally  over  the  mask. 
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a)  b) 

Figure  19.  Cross  sectional  SEM  micrographs  of  selective  are  growths  of  InP  in  stripes  oriented  along  two 
orthogonal  directions  in  the  InP  substrate. 

In  a),  the  growth  habit  shows  piling  up  at  the  edges  owing  to  rapid  surface  migration  and  limited 
(111)  plane  incorporation.  This  growth  habit  leads  to  self-limited  growth  and  is  less  useful  for 
device  applications  than  b). 

Quaternaries  behave  differently  than  InP  in  selective  area  growth.  There  is  a  stronger  tendency  for 
surface  migration  and  the  formation  of  thick  “wings”  at  the  edges  of  the  mesa.  This  is  particularly 
true  for  the  <1 10>  oriented  stripe.  Figure  20  a)  and  b)  shows  the  growth  rate  enhancement  in  the 
center  of  a  wide  stripe  for  the  both  <H0>  orientations.  The  left  picture  shows  the  growth  rate 
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Fig.  20  Growth  rate  enhancement  and  QW  emission  energy  for  MQW  layers  grown  by  SAG  on  stripes  with 
orthogonal  orientations. 


enhancement  and  the  right  shows  the  resulting  QW  emission  energy  for  a  QW  sample  selectively 
grown  in  these  openings.  The  independent  variable  in  both  cases  is  the  widths  of  the  SiNx  mask  on 
either  side  of  the  stripe.  The  wider  is  the  mask,  the  greater  is  the  growth  rate  enhancement  and  the 
longer  is  the  wavelength  of  emission.  Of  course  this  results  from  the  thicker  QW’s  grown  in  the 
stripe  opening. 
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Figure  21  shows  the  PL  spectra  of  such  samples. 
Note  that  the  linewidths  are  narrow  and  PL 
intensity  is  high.  These  data  clearly  show  that 
variable  energy  stripes  can  be  formed  by  varying 
the  width  of  the  mask  adjacent  to  the  stripes.  Such 
a  capability  allows  one  to  design  an  active  stripe 
with  engineered  energy  gaps  in  the  various  regions 
of  the  CPM  laser. 

1000  1050  1100  1150  1200  1250 

Wavelength  (nm) 

Figure  21  PL  of  QW  samples  grown  by  SAG  with 
two  different  mask  widths  showing  the  shift  to 
longer  wavelength  with  mask  width.  Note  there  is 
little  change  in  the  PL  intensity  in  SAG  grown 
materials. 


2.3  Mode  locked  laser  results 

During  the  course  of  the  program,  we  were  able  to  demonstrate  a  state  of  the  art  CPM  laser  using 
the  materials  growth  technology  described  above  to  grow  the  wafers  and  by  using  ridge  waveguide 
designs  for  the  laser.  This  was  chosen  to  minimize  the  number  of  processing  steps  in  the  early 
stages  of  the  device  studies  while  the  BH  structure  was  developed.  SAG  active  regions  presented 
an  additional  processing  difficulty  that  exacerbated  our  attempts  to  build  BH  structure  with  these 
active  regions.  The  nonplanarity  of  the  SAG  mesa,  made  it  difficult  to  define  ~lpm  wide  mesas  for 
BH  structures  on  the  mesa  top.  As  a  result  all  of  the  results  that  are  presented  in  this  section  were 
obtained  on  ridge  waveguide  structures  on  planar  substrates.  In  these  devices  an  etch  stop  layer  is 
included  in  the  structure  below  the  cap  layer.  A  2-3  micron  mesa  is  defined  and  etched  through  the 
cap  and  upper  confining  layer  to  the  etch  stop  layer  which  coincides  with  the  waveguide.  In 
addition,  small  gaps  are  etched  (10  pm  long)  in  the  cap  layer  to  provide  a  measure  of  isolation 
between  the  various  segments  of  the  CPM  device.  All  of  the  devices  described  in  this  section  were 

operated  as  passively  mode  locked  devices  as  we 
did  not  have  sufficient  RF  power  in  our 
oscillator  to  drive  the  modulator  sections  and 
had  not  yet  developed  the  strip  line  technology 
necessary  to  couple  the  ~40  GHz  RF  into  the 
modulator  sections. 

Figure  22  shows  the  L-I  characteristics  of  a 
passively  locked  CPM  laser  for  different  values 
of  the  bias  to  the  absorber.  As  the  bias  swings 
from  positive  bias  to  increasingly  greater 
negative  values,  the  threshold  of  the  laser 
increases.  This  is  to  be  expected  since  the 
absorption  at  the  laser  wavelength  in  reverse 
biased  absorber  increases  with  increasing 
reverse  bias.  The  threshold  current  increases 
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from  approximately  16  mA  at  an  absorber  bias  of  IV  to  35  mA  at  an  absorber  bias  of  -6V.  The 

quantum  efficiency  of  the  device  also  decreases  with 
increasing  reverse  bias.  The  increased  loss  with 
reverse  bias  results  from  the  quantum  confined  Stark 
effect  in  the  absorber  region.  The  absorption  of  light 
in  this  region  results  in  a  current  in  the  absorber 
contact.  This  si  illustrated  in  Figure  23,  where  the 
absorber  current  is  plotted  as  a  function  of  the  gain 
section  current.  The  presence  of  increasing  absorber 
section  current  with  reverse  bias  at  zero  gain  section 
current  is  indicative  of  leakage  between  the  sections. 


Figure  23  The  absorber  detected  current  as  a 
function  of  the  gain  section  current  in  a  CPM 
laser. 


Gradual  suppression  of  nearest  modes 
when  approaching  CPM  operation 


Devices  like  this  were  tested  under  a  variety  of  bias 
conditions.  Transient  light  characteristics  were 
measured  using  and  autocorrelator  to  measure  the 
pulse  width  under  modelocked  conditions.  The  most 
reliable  measure  of  CPM  was  the  time  averaged 
spectra.  Under  modelocked  operation,  the 


Figure  24  Spectra  of  modelocked  laser  under  three 
bias  conditions  approaching  and  entering  CPM 
operation. 

longitudinal  modes  lock  in  phase  to  produce  a 
repetitive  pulsed  output  in  the  time  domain. 

Under  symmetric  CPM  operation  the 
degeneracies  of  the  longitudinal  modes  results 
in  a  spectrum  with  only  even  modes  in  the 
spectrum,  i.e.  every  other  mode  is  missing.  A 

failsafe  signature  of  CPM  operation  was  a  spectrum  with  twice  the  mode  spacing  expected  for  the 
cavity  length.  Fig  24  shows  the  spectra  of  a  CPM  laser  under  different  bias  conditions.  Note  that 
higher  forward  bias  results  in  the  formation  of  wider  mode  spacing  characteristic  of  CPM 
operation.  Time  domain  measurements  verified  that  the  device  was  modelocked  under  those  bias 
conditions.  In  general,  the  range  of  operating  conditions  under  which  CPM  operation  was  observed 
was  narrow  for  the  devices  fabricated  in  this  program.  The  high  thresholds  observed,  the  leakage 


Figure  25  CPM  spectra  of  modelocked  laser.  The 
mode  spacing  is  twice  the  longitudinal  spacing 
expected  for  the  cavity  length.  This  is  characteristic 
of  CPM  operation. 
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between  segments  measured  and  the  use  of  ridge  waveguides  led  to  excessive  heat  dissipation.  In 
addition,  the  truncation  of  the  device  effort  owing  to  loss  of  the  student  involved  prevented  us  from 
implementing  SAG  grown  active  regions  to  “engineer”  the  cavity  for  optimal  operation.  The 
analysis  of  section  shows  that  by  balancing  the  dynamic  phase  shift  in  the  cavity  it  is  possible  to 
extend  the  range  of  passive  CPM  operation  beyond  what  has  been  previously  achieved. 

When  the  bias  conditions  are  optimized,  CPM  operation  over  a  wide  range  of  wavelengths  is 
observed.  Figure  25  a)  and  b)  shows  the  spectra  of  the  best  device  we  have  fabricated  operating 
under  the  optimal  conditions.  Double  mode  spacing  is  observed  over  a  wide  range  of  wavelengths. 
These  modes  are  expected  to  be  strongly  phase  locked  and  are  the  wavelength  sources  that  we  seek 
to  use  in  WDM  applications.  The  amplitude  of  the  modes  varies  strongly  with  wavelength.  Only 
eight  modes  have  intensities  within  a  factor  of  two  of  one  another.  Nonetheless,  these  results 
suggest  that  with  an  optimized  device  the  goals  of  the  program  can  be  realized. 


137.8mA/-5V  abs/15  deg  Cl 
2mW 


Time  (ps 

Figure  26  Autocorrelation  trace  of  CPM  laser  showing 
a  2  ps  FWHM.  This  corresponds  to  a  1.2  ps  pulse 


The  time  domain  characteristics  of  the  output 
can  be  measured  by  using  an  auto  correlation 
technique.  Here  the  pulse  is  split  and  counter 
propagated  through  a  doubling  crystal.  The  delay 
between  the  arrival  of  the  two  halves  of  the  pulse 
is  varied  to  allow  them  to  interact  in  a  region  of 
the  crystal  under  a  slit.  By  varying  the  delay  and 
measuring  the  second  harmonic  radiation 
through  the  slit  as  a  function  of  the  delay  one  can 
measure  the  correlation  of  the  pulse  with  itself. 
The  resultant  autocorrelation  output  shown  in 
Figure  26  shows  a  pulse  width  of  2  ps. 
Deconvolving  the  pulses  from  one  another,  we 
calculate  a  pulse  width  of  1 .2  ps.  This  is  equal  to 
the  shortest  CPM  pulse  observed  in  the 
literature  and  is  consistent  with  the  measure 
linewidth  of  the  laser  mode. 
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3.0  Theory  of  Passively  Mode-Locked  Semiconductor  Laser  under  Self  Phase 
Modulation 


3.1  Introduction 

There  has  been  growing  interest  in  the  passive  mode  locking  of  monolithic  semiconductor  lasers  to  produce 
higher  repetition  frequency  light  sources  than  allowed  by  the  electrical  modulation.  High  repetition 
frequencies  ranging  from  hundreds  of  GHz  up  to  more  than  THz  with  sub-picosecond  pulses  have  been 
demonstrated  by  shortening  the  cavity  length  [1],  by  employing  multiple  colliding  pulse  mode-locking  [2], 
or  by  employing  harmonic  mode-locking  [3].  Compact  in  size  and  self-aligned,  the  passively  mode-locked 
laser  diode  operating  at  20-100  GHz  frequency  has  potential  of  the  light  source  without  electrical  parasitic 
limitations  in  the  wavelength  division  multiplexing  (WDM)  system.  Also,  the  multiple  wavelength 
components  in  the  mode-locked  pulse  are  inherently  phase-coherent  and  equally  spaced  in  the  spectrum, 
which  can  open  a  new  application  in  WDM  system  design  [4],  [5]. 

To  understand  the  mode-locking  mechanism  where  the  absorption  and  gain  saturation  by  an  optical 
pulse  and  their  recovery  interplay  with  each  other,  device  simulation  is  indispensable.  In  many  cases,  the 
electric  field  propagation  through  the  multi-section  cavity  is  simulated  in  time  domain  with  the  spontaneous 
emission  as  a  seed,  then  the  steady  state  after  many  round  trips  is  investigated  [6]-[l  1].  However,  for  the 
monolithic  device  design  with  a  high  repetition  frequency  (>50  GHz),  the  mode-locked  spectrum  contains 
only  several  modes  separated  by  more  than  a  few  A’s.  Also,  at  high  frequency,  the  complete  recovery  of  the 
gain  or  absorption  saturation  is  not  allowed,  which  means  that  the  carrier  density  modulation  is  not  as  deep 
as  the  optical  modulation.  In  this  case,  a  frequency  domain  simulation  using  the  multiple  modes  as  the  basis 
is  more  attractive  than  the  time  domain  simulation  that  can  include  the  spectral  dependencies,  at  best,  by 
gain-spectrum  filtering  and  the  group  velocity  dispersion  [7],  [9].  The  frequency  domain  method  has  been 
used  either  for  qualitative  descriptions  [2],  [12]  or  in  a  semi-quantitative  way  including  only  the  nearest 
neighbor  mode-coupling  [13].  In  this  work,  passively  mode-locked  laser  diodes  were  simulated  using  the 
frequency  domain  method  that  includes  higher  order  mode  coupling  than  the  nearest  neighbor  coupling  and 
the  carrier  density  dependencies  of  the  gain,  carrier  recombination  rate,  and  refractive  index  change.  Our 
recent  study  [14]  showed  that  the  specific  carrier  density  dependencies  of  the  material  gain  and  the  carrier 
recombination  rate  are  the  key  parameters  determining  the  gain  saturation.  The  line  width  enhancement 
factor  is  also  known  to  be  a  function  of  the  carrier  density  and  the  operation  conditions  [15].  Therefore,  the 
dependencies  of  the  differential  gain,  carrier  lifetime,  and  the  line  width  enhancement  factor  on  the  amplifier 
current  or  absorber  reverse  bias  should  be  taken  into  account  to  predict  the  device  characteristics  over  a  wide 
operation  range.  We  will  show  how  the  inclusion  of  the  self  phase  modulation  (SPM)  accompanying  the 
mode  coupling  reduces  the  possible  operation  range  of  the  mode  locking  and  also  why  the  most  of  the 
experimental  results  to  date  show  transform  limited  pulse  width  even  under  SPM.  The  rate  equation  for  the 
supermode  spectrum  is  derived  in  the  section  2  with  the  emphasis  on  the  carrier  density  dependence  of  the 
parameters.  The  section  3  explains  the  self-consistency  of  the  solution,  numerical  procedure  for  handling  the 
SPM,  and  the  solution  stability.  The  device  characteristics  simulated  under  various  operation  conditions  are 
presented  in  section  4,  which  is  followed  by  section  5  for  summary. 


3.2  Rate  equations  for  a  mode-locked  laser  diode 

Schematically  shown  in  Fig.  27  is  the  model  device  with  a  cavity  extending  over  the  range  0 <z<L.  The 
absorber  section  of  length  haL  is  centered  at  (1  +  0)/.  /  2  .  The  rest  of  the  cavity  is  assumed  to  be  amplifier 
section.  Here,  we  assume  that  the  total  absorption  inside  the  absorber  is  not  high  and  both  facets  of  the 
cavity  are  either  as  cleaved  or  high-reflection  coated.  Under  these  conditions,  the  average  photon  density  can 
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be  approximated  to  be  constant  over  the  whole  cavity  when  the  intensities  of  forward  and  backward 
propagating  waves  are  summed  [16]. 


Fig.  27.  Schematic  figure  showing  the  device  dimensions. 


The  electric  field  is  given  by 


E(jr,t)  =  Re  ^ Em  (t)e,mJa(x, y)Jj  cos(p, 


where  a(x,y)  represents  the  cross-sectional  field  distribution  and  the  wave  vector  (3m  satisfies  the  relation 
of  Pm+I  ~(im  =n/L.  Here,  the  w-th  mode  electric  field  Em  is  a  complex  number  and  slowly  varying  in 
time.  When  the  device  is  uniformly  pumped  at  its  threshold,  the  optical  angular  frequency  is  given  as  2itvm , 
which  is  the  Fabry-Perot  mode  of  the  cavity.  When  the  device  is  in  mode-locked  operation,  the  optical 
frequencies  may  be  shifted  from  their  original  values  but  they  are  still  expected  equally  spaced.  When  80  is 
the  shift  of  the  central  mode  (tn= 0)  that  we  arbitrarily  assigned  at  the  gain  peak,  the  shift  of  the  other  modes 
can  be  described  by  co„,  =  2nvm  +  50  +  mb  so  that  the  repetition  frequency  Q  is  given  by  the  relation  of 


Q  =  co„ 


-vm)  +  8. 


The  temporal  and  spatial  behavior  of  the  photon  density  in  the  active  gain  or  absorber  can  be 
obtained  by  taking  the  average  of  \E(r,t)\2  spatially  over  a  few  wavelengths  and  temporally  over  several 
optical  oscillations.  This  spatial  average  implies  the  smearing  of  the  index  grating  due  to  carrier  diffusion. 
From  (1),  the  z  and  t  dependence  of  the  photon  density  S  is  given  by 


S(z,t)  =  S0  1  +  X2 


!cos(  — |Re[s,(*)eto' 


where 


C k )  _  Em+\Em 


I  mBmK 

represents  the  depth  of  the  optical  field  modulation  at  frequency  kQ.  To  clarify  the  meaning  of  S0,  the 
average  cross-sectional  optical  power  is  given  as  EphvgS0(wd)/T ,  where  Eph  is  the  photon  energy,  vg  is 

the  group  velocity  of  the  light,  (wd)  is  the  cross-sectional  area  of  the  active  region,  and  T  is  the  optical 
confinement  factor. 
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The  photon  density  modulation  at  various  frequencies  kQ  (k  =  1,2, . . .)  induces  a  carrier  density 
modulation  at  those  frequencies  described  as 


n  =  n  +  ^  2  cosl 
*>o  V 


(5) 


Throughout  the  paper,  a  quantity  with  a  bar  means  the  average  over  the  time  (or  DC  quantity)  and  that  with 
a  tilde  such  as  n(k)  or  7(k)  represents  the  modulation  amplitude  at  the  frequency  kQ.  It  should  be  noticed 
that  these  quantities  are  different  in  the  absorber  and  in  the  amplifier  regions. 

In  the  amplifying  section  under  current  pumping,  the  carrier  density  is  related  to  the  photon  density 
through  the  following  rate  equation 


R(n)-v  Sg(n) 


(6) 


dn_  _  J non 

dt  q 

where  J nom  is  the  volume  current  density,  q  is  the  electron  charge,  R{n)  is  the  carrier  recombination  rate, 

and  g(ti)  is  the  optical  gain.  If  the  carrier  modulation  terms  « (i)  are  small  enough,  R{n)  and  g(n)  can  be 
linearly  expanded  around  n  =  n  as  in 


and 


Rjn)  =  Rjn )  +  —  ^ 2 cos  — -  Re[n  (ArVin/  ] 
TkZo  L  / 

k>  0  \  L  J 


(7) 


(8) 


where  1  /  x  =  (dR  /  dn)  ^  and  G  =  ( dg  /  dn)  | ^ .  It  should  not  be  overlooked  that  x  and  G  strongly  depend  on  n 
which  is  also  changed  by  S0 . 

The  substitution  of  (3),  (7),  and  (8)  in  (6)  and  reorganizing  the  terms  up  to  the  first  order  of  n(k)  and 


S' W  gives  the  DC  operation  point  determined  by 


J „ 


-R(n)-v  gS0  =0 


q 

and  the  relationship  between  the  carrier  and  photon  density  modulations  as 


-gS  0 


r(k) 


(9) 


(10) 


(GS0  +\/x)  +  ikQ 

For  the  absorber  section,  the  photo-generated  carriers  are  swept  away  by  the  reverse  bias,  which  is  better 
described  by  the  phenomenological  time  constant,  i„ .  The  introduction  of  t„  changes  (9)  to 


- +  R(n)  =  v  i-g)S0 

and  in  (10),  the  carrier  lifetime  t  in  the  absorber  section  should  be  redefined  as 
1 


1  (a  bs) 


1  dR 
x,u,  dn 


(ID 


(12) 


These  changes  make  the  following  derivations  applicable  to  the  absorber  section  as  well  as  the  gain  section. 
In  the  high  repetition  rate  cases  that  we  are  mostly  interested  in,  Q  is  so  much  larger  than  GS0  or  l/x  that  the 
carrier  modulation  usually  lags  the  optical  modulation  by  almost  %!2  and  is  also  small  in  magnitude,  which 
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validates  the  linear  expansion  in  (7)  and  (8).  The  carrier  recovery  in  the  reverse  bias  in  the  saturable 
absorber  is  described  based  on  the  experimental  results  [17],  [18]. 

The  gain  and  the  refractive  index  modulation  due  to  the  carrier  density  modulation  given  in  (5)  is 
the  origin  of  the  mode  coupling  between  the  neighborhood  modes  and  the  self  phase  modulation, 
respectively.  The  effect  of  modulating  optical  gain  and  the  refractive  index  on  a  mode  can  be  obtained  by 
the  overlap  integral  of  the  mode  with  the  polarization  [13].  This  overlap  integral  can  be  separated  into  the 
cross-sectional  integral  and  the  subsequent  longitudinal  one  along  the  cavity  for  heuristic  purposes.  The 
cross-sectional  overlap  integral  resulted  in 


X  T,~i('2nVm  _®m)+ rT 


cos(p  mz)e^  =£rvf 


gi(z,0  .©/ 


i—An,(z,t)  El  cos(P,z)e  '  , 
c 


where  xpm  is  the  photon  lifetime  at  co„, ,  which  is  determined  by  the  internal  loss  and  the  distributed  mirror 
loss,  and  Ap.,  represents  the  refractive  index  change  at  co,  from  its  laser  threshold  value.  Also,  the  spectral 
dependence  of  the  gain  function  is  implemented  by  g{n)  =  gt  (n) .  The  optical  confinement  factor  T 


reflects  the  fact  that  the  gain  and  the  refractive  index  change  only  in  the  active  region  of  the  cross-section. 
Using  the  definition  of  the  line  width  enhancement  factor 


Jdyjdv)  (14) 

(di]ldn). 

where  X/  =  X/  +  ‘li  is  the  electric  field  susceptibility  at  ©, ,  the  gain  and  refractive  index  change  for  the  /-th 
mode  can  be  described  as 


2  c  2  J  2 

where  g,  =g,(n)  and  ic;  =  (o/Ap/(«)/c .  In  (15),  most  of  the  quantities  are  position  dependent  implicitly. 

Substituting  (5),  (10),  and  (15)  in  (13)  then  reorganizing  the  terms  according  to  their  angular 
frequencies  makes  it  possible  to  find  the  rate  equation  for  each  mode  Em  which  is  coupled  to  the 
neighboring  modes  given  by 


=  -  z'(80  +  mb)  ■ 


where  the  mode  coupling  efficiency  C ^  and  the  modulation  field  distribution  D^k\z)  defined  as 


(GS0  + 1  /  x)  +  ik£l 


D(k)(z)  =  -  +  -  cos 
2  2 


(18) 


respectively.  It  is  noticed  that  the  definitions  of  s(k)  and  C(„k)  are  extended  to  negative  k  indices  too  by 
using  s(~k)  =  [s(A)]‘  and  C(mk)  =  [C^0]’ .  The  spatial  average  of  (16)  over  the  cavity  length  L  completes  the 
derivation  of  the  rate  equation  for  Em ’s.  When  the  device  cavity  is  divided  into  the  absorber  and  amplifier 
sections  whose  characteristics  are  constant  in  each  section,  the  spatial  average  can  be  done  very  easily  by 
multiplying  ha ,  the  fraction  of  the  cavity  occupied  by  the  absorber,  to  the  absorber  parameters  and  (\-ha )  to 
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the  amplifier  parameters.  The  explicit  position  dependence  of  Di"k\z)  can  be  reflected  by  the  modification 
of 

=f 

and 

<20) 

where  0  is  the  parameter  showing  the  deviation  of  the  absorber  section  from  the  center  of  the  cavity  as 
depicted  in  Fig.  1.  The  rate  equation  for  the  multimode  system  can  be  expressed  in  the  matrix  form  as 


knh„ 


(19) 


dE» 

dt 


ml  ^1  > 


(21) 


i 

where  Aml  is  the  corresponding  term  in  (16)  after  averaging  overZ.  Simultaneous  lasing  of  many  coupled 
modes  can  be  a  steady  state  solution  of  system  when  the  matrix  A  in  the  rate  equation  (21)  has  the 
eigenvalue  of  zero,  which  is  called  the  supermode.  In  the  simplest  case,  k  =  ±  1  coupling  is  included,  the 
absorber  is  located  at  one  end  of  the  cavity  by  setting  0  =  1  -  ha ,  and  g  =  G(n  -n0)  and  R  =  n/x  for  all 
carrier  density  ( n0  is  the  transparency  carrier  density),  then  the  steady  state  of  the  rate  equation  (21)  can  be 
reduced  to  the  Eq.  (33)  of  Ref.  [13]. 


3.3  Simulation  procedure 

Finding  the  steady  state  solution  of  (21)  is  restated  as  an  eigenvector  problem  to  find  a  zero  eigenvalue,  but 
solving  the  problem  is  much  more  complicated  than  it  apparently  looks.  First  of  all,  the  set 

{S'(i)}(Z  =  l,2,...)  obtained  from  the  solution  {£,„ }  using  (4)  should  be  matched  to  the  corresponding 
values  used  in  A.  This  self-consistency  of  {7(k)}  is  the  major  difference  of  this  problem  from  the  ordinary 
eigensystems.  Among  the  various  eigen-solutions  obtained  from  a  given  A,  only  the  one  giving  self- 
consistent  is  meaningful.  In  numerical  pursuit  of  the  solution,  we  keep  following  the  eigenvector  that 

has  the  maximum  real  part  of  the  eigenvalue,  and  then  try  to  adjust  parameters  until  this  eigenvalue  becomes 
zero  in  both  real  and  imaginary  parts.  The  choice  of  the  maximum  real  part  of  the  eigenvalue  is  proposed  on 
the  basis  that  the  eigenvalue  of  A  is  interpreted  as  the  change  rate  of  the  supermode  according  to  (21),  so  the 
supermode  with  the  maximum  real  part  of  eigenvalue  will  dominate  finally. 

For  the  self-consistency  of  {7^} ,  the  phase  of  T<i!  is  as  important  as  its  magnitude.  It  seems  that  the  phase 
of  one  7l'k'>  can  be  arbitrarily  changed  by  the  choice  of  another  time  reference  which  changes  the  relative 
phase  between  Em 's.  However,  the  use  of  a  different  time  reference  changes  the  phases  of  7(k)  in  (3),  so  the 
solution  {Em }  cannot  generate  the  same  phase  of  7{k)  as  the  input.  Therefore,  the  self-consistency  of  the 
7^  phase  can  be  achieved  only  by  finding  the  right  solution  parameters.  Still,  to  eliminate  the  ambiguity  in 
choosing  the  time  reference,  the  phase  of  one  7 (*)  is  set  zero.  We  usually  choose  kn,f  =  1  for  this  purpose, 
but  the  choice  of  kref  is  not  very  crucial  when  several  k  numbers  are  included  in  the  simulation. 

In  the  simulation  procedure,  the  average  photon  density  S0  is  the  only  independent  variable.  All  the 
other  parameters  are  interrelated,  so  they  should  be  solved  self-consistently  by  iteration  as  well  as  numerical 
solution  searching  procedures.  Each  iteration  step  starts  with  a  set  of  Em  which  is  given  either  by  a  guess  at 


Monolithic  Modelocked  Lasers  for  WDM 
Dapkus/USC 


18 


May,  2001 


first  or  as  the  eigenvector  solution  of  A  from  the  previous  iteration  step.  The  correlation  terms  {s for  A 
may  be  preset  using  this  {Em} .  However,  at  the  very  beginning  of  the  iteration  procedure,  the  parameters  do 
not  generate  a  self-consistent  set  of  {. s(k) }  very  well.  In  this  case,  it  is  recommended  to  handle  these  s(k)'s 
as  independent  parameters  instead  of  as  functions  of  Em' s.  That  is,  {v  ■*'  ' }  is  arbitrarily  adjusted  like  the 
other  parameters  to  minimize  the  difference  between  the  given  (s(k))  and  the  correlation  calculated  by  (4) 
using  the  solution  eigenvector  for  Em  values. 

When  {£„,}  and  {?(<r)}  are  given,  the  average  carrier  density  in  the  absorber  section,  n(ahs) ,  is 
obtained  by  (1 1),  in  which  the  average  absorption  [-g(aAv)]  depends  not  only  on  n(ahs)  but  also  on  Em  since 
the  average  gain  g  should  be  defined  as 


g  = 


U  E- 1! 


(22) 


to  reflect  the  spectral  dispersion  of  gm  in  the  carrier  density  rate  equation  which  was  oversimplified  in  (16). 
The  average  differential  gain  G,  which  was  used  in  (8),  (10),  and  (17)  without  clear  definition,  should  be 
defined  in  a  similar  way.  Once  n^abs^  is  found  for  the  given  { Em } ,  all  the  parameters  in  the  absorber  section 


can  be  determined  from  their  definitions.  The  average  carrier  density  in  the  amplifier  section,  n^amp^ ,  cannot 


be  directly  obtained  in  the  same  way  as  n{ahs)  since  Jnom  is  unknown  yet.  Instead,  n{amp)  should  be  guessed 

to  set  all  the  gain  section  parameters,  then  the  80  and  8  parameters  should  be  picked  up  to  determine  the 
matrix  A  completely.  The  eigenvalue  with  the  maximum  real  part  and  the  corresponding  eigenvectors  are 
obtained  from  A,  and  the  resulting  {s(k)}  is  also  determined.  By  iterating  n{amp) ,  80,  and  8  parameters, 

these  three  parameters  can  be  best  fitted  to  make  both  the  eigenvalue  and  the  phase  of  s  ”{  as  small  as 
possible.  The  phase  of  s{k'c/)  is  mostly  affected  by  the  variation  of  the  5  parameter,  while  n(amp)  and  80 
mostly  affect  the  real  and  imaginary  parts  of  the  eigenvalue,  respectively.  When  these  three  parameters  are 
best  adjusted  to  minimize  the  eigenvalue  and  the  phase  of  the  5  ,  the  new  set  of  {Em }  is  obtained  from 

the  eigenvector,  which  completes  one  iteration  step.  The  iteration  process  should  be  repeated  until  the  given 
{Em}  reproduces  itself.  When  the  self-consistent  solution  is  obtained  for  the  given  S0,  the  amplifier 
injection  current  required  to  maintain  the  supermode  can  be  easily  calculated  from  the  Jnnm  in  (9).  Also  the 
pulse  width  can  be  deduced  by  the  Fourier  transform  of  {•Z^e^0’"-010^ }  to  the  time  domain.  By  continuously 
varying  S0  value  then  obtaining  the  solution  including  Jnnm ,  the  device  performances  can  be  predicted  as 
functions  of  the  amplifier  current. 

It  should  be  remembered  that  there  always  exists  an  extra  trivial  solution  that  is  the  single  mode 
lasing  whose  threshold  gain  for  the  amplifier  section  is  the  sum  of  the  cavity  loss  and  the  absorption  in  the 
absorber  at  the  lasing  wavelength.  When  a  supermode  solution  is  found  with  the  amplifier  carrier  density 
n(amp) ,  it  can  be  an  actual  solution  only  when  the  gain  of  any  single  mode  at  n(amp)  is  not  higher  than  its 

threshold.  Therefore,  the  magnitude  of  the  threshold  modal  gain  reduction  with  respect  to  the  single  mode 
lasing,  which  is  defined  as 


Rg  -  mini  +  [  r Sm(ahs)^a  ^Sm(amp)^  >  (^3) 

IAV  J 

is  a  good  indicator  of  how  much  the  supermode  is  more  stable  than  any  possible  single  mode  lasing  in  the 
system. 
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In  a  preliminary  study,  it  was  found  that  a  self-consistent  solution  is  very  hard  to  find  at  some  S0 
values  when  the  line  width  enhancement  factor  is  not  negligible.  In  contrast,  the  simulation  of  an  imaginary 
device  that  does  not  have  any  refractive  index  change  shows  stable  solution  over  a  wide  range  of  S0  value. 
The  dotted  line  in  Fig.  28  shows  one  of  the  typical  examples.  As  the  photon  density  increases,  first,  the 
absorption  saturation  grows  to  increase  the  stability  of  the  mode-locked  pulse  as  can  be  seen  in  Fig.  2-(a). 
As  the  photon  density  increases  further,  the  increase  of  the  absorber  carrier  density  reduces  the  differential 
absorption.  As  a  result,  the  contribution  of  the  gain  saturation  starts  to  compete  with  that  of  the  absorption 
saturation,  and  beyond  a  certain  point,  further  increase  of  the  photon  density  starts  to  decrease  the  mode¬ 
locking  stability.  Also,  because  of  the  time  lag  between  the  photon  density  and  the  carrier  density  as  is 
quantitatively  described  in  (10),  the  absorption  and  gain  saturation  introduces  phase  shift,  which  can  result 
in  the  change  of  the  repetition  frequency  (5)  even  without  the  SPM  effect.  Of  course,  in  the  absence  of  the 
SPM,  the  chirp  parameter  defined  as  the  second  time  derivative  of  the  phase  of  the  electric  field  in  time 

domain  is  zero.  The  dotted  lines  in 
Fig.  28-(b)  and  -(c)  show  these 
0.5  r - -  |  trends.  Therefore,  the  passive  mode- 

(a)  locking  without  the  SPM  effect  is 

0.4-  ,  '  *  \  described  by  the  competition 

^  '  '  between  the  absorption  and  the  gain 

'g  0.3  -  /  \  saturation,  where  the  carrier  density 

3  '  \  response  in  phase  with  the  photon 

of  0.2  ~  *  density  modulation  determines  the 

g  \  stability  of  the  mode-locking,  while 

0.1  “  if  *  the  out  of  phase  contribution 

,  .  ,  determines  the  repetition  frequency 

■u  ’  tuning.  Even  though  not  exactly 

(b)  200  -  /  matched,  the  photon  density  where 

^  Ann  '  the  mode-locking  stability  or  R„ 

n  TUU  -  , 

^  ,  becomes  maximum  is  close  to  that  of 

0  “  /  zero  8,  as  they  are  determined  by  the 

cn  "m-  , '  balance  of  real  and  imaginary  part, 

-1 00  -  '  respectively,  between  the  absorber 

and  amplifier  contributions. 


-  -  -  y  =  0 

- y  allowed 


/_\  1.0“  _  _  _  _  q  Fig.  28.  Effect  of  the  line  width 

'  '  Y  X0  8  _  ’  ll  H  enhancement  factor  on  the  parameters  of 

fjP1  q  5  _  abs  ^  3  0We  (a)  mode-locking  stability  Rg ,  (b) 

3  /  repetition  frequency  shift  8/2 n,  and  (c) 

^  0.0  —  y  - 8 -  chirp  parameter  d2§/dt2  as  functions 

.  JP  of  the  photon  density.  Dotted  and  the 

cn®"  _0.5  -  x-|  2  solid  lines  are  simulated,  respectively,  in 

^abs  ’  the  absence  and  presence  of  the  line 

.  ,  ,  ,  |  |  ,  i  ,  ,  ,  |  width  enhancement  factor,  while  the 

•  o  4  16  36  64  markers  shows  the  cases  where  T  in  the 

absorber  section  is  arbitrarily  increased 
^  ^  A  or  decreased  from  the  actual  value.  The 

Optical  Power  Flux  in  the  Cavity  (mW)  optical  power  is  plotted  in  a  Square.root 

scale  to  clarify  the  features  all  over  the  range.  The  simulation  is  performed  for  a  strained  4  quantumwell  active  region 
with  the  parameters  of  L=0.5  mm  (Q~80  GHz),  ha  =0. 1  (located  at  one  end),  and  ^20  ps. 


Yabsx1-2 
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The  existence  of  the  solution  over  a  wide  range  of  S0  when  the  SPM  is  not  allowed  can  be  explained  from  a 
mathematical  point  of  view  as  follows.  Keeping  it  in  mind  that  Gm  depends  on  m  very  weakly,  the  matrix^ 
is  almost  Hermitian  matrix  when  Km  =  0  and  yOT  =  0 .  Therefore,  the  eigenvalues  are  real  and  can  be  sorted 
according  to  their  magnitude,  which  makes  it  unambiguous  to  pick  up  the  solution.  When  the  SPM  is 
allowed,  its  contribution  can  be  considered  in  two  ways.  When  the  contribution  due  to  the  dynamic 
modulation  of  the  carrier  density  is  neglected  (y„,  =0),  the  inclusion  of  i[ha  Km(afc)  +0  ~K)Km(amP)  'n  the 

diagonal  elements  of  A  can  be  just  cancelled  by  a  small  change  of  80  and  8,  which  makes  A  Hermitian-like 
again.  This  is  possible  as  the  refractive  index  change  as  a  function  of  the  frequency  can  be  approximated  as 
a  linear  function  when  the  mode-locked  spectrum  is  much  narrower  than  the  gain  spectrum  itself.  The  small 
change  of  80  and  8  is  the  Fabry-Perot  mode  shift  from  the  uniformly  pumped  laser  as  a  result  of  the  net 
index  change  after  a  round  trip  through  the  amplifier  and  the  absorber  sections. 

However,  when  ym  is  allowed,  the  matrix  A  is  no  longer  Hermitian-like  and  there  may  exist  a 

couple  of  eigenvalues  whose  real  parts  are  close  to  each  other.  In  this  case,  because  of  the  supermode 
competition,  a  small  change  in  the  device  parameter  can  drive  the  simulation  into  a  completely  different 
result.  When  the  dynamic  SPM  is  allowed  through  nonzero  y„, ,  a  self-consistent  and  stable  solution  can  be 
obtained  only  when  the  contribution  from  the  absorber  section  is  partially  cancelled  by  that  of  the  amplifier 
section.  This  is  possible  when  the  term  gy „,C^ Z)(i)  of  the  absorber  and  the  amplifier  sections  are  almost 
opposite  to  each  other.  In  the  case  of  the  SPM-free  case,  the  balance  between  gC^D^  terms  allowed  the 
maximum  Rg  and  zero  8.  As  the  line  width  enhancement  factor  for  the  absorber  is  smaller  than  that  of  the 

amplifier  because  of  the  difference  in  the  differential  absorption  and  gain,  the  balance  between  the  dynamic 
SPM  can  be  achieved  at  lower  power  than  the  maximum  Rg  condition  of  the  SPM-free  case.  When  the 

dynamic  SPM  is  not  exactly  balanced,  the  solution  stability  is  at  risk  and,  beyond  a  certain  point,  the 
solution  changes  discontinuously  as  the  optical  power  changes  continuously,  which  is  the  point  where  the 
SPM  solution  (solid  lines)  is  not  drawn  any  further  in  Fig.  2.  In  the  figure,  the  Km  parameters  are  set 
arbitrarily  zero  to  facilitate  the  comparison.  Except  the  chirp  characteristics  in  Fig.  2-(c),  the  Rg  and  8 

parameters  are  not  changed  much  from  the  SPM-free  solution,  because  the  solution  is  obtained  after  the 
dynamic  SPM  is  balanced  between  two  sections.  Also,  the  chirp  characteristics  are  changed  from  down- 
chirp  to  up-chirp  as  the  optical  power  increases.  When  the  absorber  line  width  enhancement  factor  is 
arbitrarily  decreased  by  20  %,  it  is  clearly  seen  that  the  reducing  the  dynamic  SPM  effect  of  the  absorber 
section  shifts  the  operation  range  to  the  lower  power  side  and  the  pulse  becomes  more  up-chirped.  The 
increase  of  the  absorber  dynamic  SPM  induces  the  opposite  result.  The  inclusion  of  the  dynamic  SPM  in  the 
passively  mode-locked  laser  considerably  narrows  the  operating  range  to  the  lower  power  side  than  the  most 
stabilized,  SPM-free  case,  and  the  stronger  absorber  dynamic  SPM  allows  more  down-chirp  and  higher 
power  operation. 

In  the  simulation  procedure,  the  carrier  lifetime  in  the  absorber  section  is  mostly  determined  by  x  w , 

which  is  not  only  the  strong  function  of  the  reverse  bias  applied  on  the  absorber  section  but  also  influenced 
by  the  photo-generated  carrier  screening.  Therefore,  the  simulation  over  a  wide  range  of  optical  power  with 
fixed  xw  is  not  completely  accurate  for  the  case  of  a  constant  absorber  reverse-bias.  Therefore,  the 
simulation  result  plotted  as  a  function  of  the  amplifier  current  with  fixed  t„  may  not  represent  the  device 
result  for  a  constant  absorber  bias.  However,  as  the  operation  range  becomes  very  limited  as  a  result  of  the 
dynamic  SPM,  using  a  constant  xiW  for  a  fixed  absorber  bias  can  be  justified  in  the  following  simulation. 
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3.4  Simulation  Result 


Using  the  simulation  method  described  above,  we  calculated  the  device  characteristics  of  passively  mode- 
locked  lasers  of  various  device  structures  under  various  operating  conditions.  As  mentioned  in  the  previous 
section,  there  exists  a  limited  operation  range  in  which  the  simulation  gives  a  stable  solution.  For  example, 
the  absorber  recovery  time  strongly  affects  the  range  of  the  laser  current  for  a  stable  mode  locking.  Figure 
29  shows  the  simulation  result  of  the  average  laser  output  power  as  a  function  of  the  injection  current  for  a 
laser  with  a  cavity  length  of  0.5  mm,  where  10%  of  laser  cavity  is  used  for  a  saturable  absorber.  By 
decreasing  the  sweep  time  constant  or  by  applying  more  reverse  bias  on  the  absorber  section  in  the  real 
device  case,  the  photo-generated  carrier  density  in  the  absorber  section  becomes  smaller  when  a  certain 
amount  of  light  power  exists.  This  means  that  absorption  saturation  becomes  comparable  to  the  gain 
saturation  at  higher  optical  power  than  the  longer  sweep  time  case.  As  the  stable  mode  locking  is  allowed 
only  when  the  chirp  due  to  absorber  section  balances  that  of  the  gain  section,  the  shorter  the  sweep  time  is, 
the  higher  is  the  allowed  optical  power.  Of  course,  high  mode-locked  laser  power  requires  high  injection 
current.  Therefore,  in  order  to  produce  a  mode-locked  pulse  power  that  is  high  enough  for  real  applications, 
it  is  required  to  apply  large  reverse  bias  to  accelerate  sweeping  the  photo-generated  carriers  out  of  the 
saturable  absorber. 


Current  (mA) 

Fig.  29.  Simulated  mode-locked  optical  power  vs.  current  relations  for  different  absorber  sweep  time 
constant. 

The  result  in  Fig.  29  can  be  rearranged  to  show  the  operation  current  range  as  a  function  of  inverse 
sweep  time  (l/tsw)  that  is  given  in  Fig.  30.  Here,  1/tsw  is  roughly  proportional  to  reverse  bias  on  the  absorber 
section.  This  simulation  shows  almost  the  same  features  as  the  experimental  results  in  Ref.  [1].  As  can  be 
seen  in  Fig.  28-(a),  the  low  current  side  of  the  stable  mode-locking  operation  current  is  mostly  limited  by 
more  stable  single  mode  lasing  (i?g<0).  Therefore,  as  the  injection  current  increases,  the  device  operation 
scheme  will  change  from  a  single  mode  lasing  to  a  multi-mode  one  which  is  the  mode-locked  pulse. 
Meanwhile,  the  high  current  side  of  the  mode-locked  solution  is  not  limited  by  the  competition  with  any 
single  mode  lasing.  In  the  simulation  with  SPM  for  the  high  optical  power,  there  exist  competing  multi- 
mode  solutions  that  are  different  in  amplitude  and  phase  but  are  very  close  in  threshold  gain.  Therefore,  any 
small  change  in  the  real  device  operation  can  shift  the  operation  point  from  one  solution  to  another.  Since  all 
of  them  are  mode-locked  pulses,  it  is  expected  that  pulses  will  be  generated  from  time  to  time  but  without 
any  clear  correlation  to  each  other.  This  is  the  Q-switched  pulse  generation.  This  explanation  well  matches 
to  the  experimental  result  given  in  Ref.  [1]. 
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Fig.  30.  Simulated  device  operation  conditions  giving  stable  mode-locked  pulse. 

Because  the  simulation  procedure  produces  the  multi  mode  spectrum  of  the  mode-locked  pulse,  the 
pulse  width  can  be  calculated  by  a  simple  Fourier  transformation  of  the  given  spectrum.  Figure  31  shows  the 
minimum  pulse  width  that  is  obtained  by  adjusting  the  injection  current  at  each  sweep  time  constant.  In  the 
allowed  operation  range  of  the  injection  current  for  a  given  sweep  time  constant,  the  pulse  width  depends 
very  weakly  on  the  injection  current  for  most  of  the  cases  investigated  in  this  study.  This  is  partly  because, 
within  the  allowed  injection  current  range,  the  variation  of  the  absorber  carrier  density  is  not  large  enough  to 
change  the  relative  importance  between  the  absorber  and  gain  sections.  However,  by  decreasing  the  sweep 
time  constant,  the  operation  optical  power  increases  and  the  absorber  saturation  becomes  more  important  in 
determining  the  operational  balance  point.  For  a  given  device  geometry,  narrower  pulses  can  be  generated 
by  decreasing  the  sweep  time  constant.  In  a  very  simplified  time-domain  explanation  of  the  mode-locked 
pulse,  the  role  of  the  saturable  absorber  is  described  as  absorbing  the  very  beginning  part  of  the  pulse  and 
being  saturated  then  letting  the  rest  of  the  pulse  penetrate  it  after  saturation,  which  cuts  off  the  rising  edge  of 
the  pulse.  In  the  frequency-domain  description,  the  threshold  gain  reduction  by  correlated  multi-mode 
locking  is  the  result  of  absorption  reduction  in  the  saturable  absorber.  The  more  number  of  modes  are  locked 
to  each  other,  the  narrower  becomes  the  pulse.  In  order  to  generate  narrower  pulse,  it  is  required  to  drive  the 
device  into  the  condition  that  allows  more  absorption  saturation.  Both  descriptions  emphasize  the  role  of  the 
strong  absorption  saturation.  Extending  this  idea  further,  we  compared  two  device  structures  with  different 
absorber  lengths.  In  Fig.  31,  two  mode-locked  devices  of  0.5  mm  cavity  length  are  simulated  where  the 
absorber  lengths  are  5  and  1 0  %  of  the  cavity  length,  respectively.  It  is  observed  that  the  device  with  the 
longer  absorber  section  generates  narrower  pulses  than  the  shorter  absorber  length  device  does.  For  an 
extended  cavity  mode-locked  laser  operated  at  a  low  repetition  frequency,  it  is  believed  that  the  absorber 
length  should  be  short  to  generate  a  narrow  pulse.  This  is  based  on  the  postulation  that  the  well-controlled 
narrow  pulse  will  be  spatially  confined  within  the  absorber  section  length.  This  might  be  true  when  tens  of 
picoseconds  width  pulses  are  generated  from  a  device  with  long  enough  absorber.  However,  for  the 
passively  mode-locked  device  operated  at  a  very  high  frequency  of  ~100  GHz,  the  pulse  period  is  already 
comparable  to  the  pulse  width  so  that  the  realistic  pulse  is  hardly  confined  inside  of  the  very  short  absorber. 
In  this  case,  the  minimum  pulse  width  is  mostly  limited  by  the  amount  of  attainable  absorption  saturation, 
which  explains  why  the  longer  absorber  is  more  suitable  to  generate  narrow  pulses. 
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Fig.  3 1 .  The  minimum  pulse  width  attainable  with  the  given  sweep  time  constant  for  two  different  absorber 
length  cases. 


3.5  Conclusion 

We  developed  a  frequency-domain  simulation  method  that  is  suitable  for  a  high  frequency  passive  mode 
locking.  This  method  includes  the  specified  carrier  density  dependencies  of  the  gain  and  absorption  spectra, 
and  also  the  self  phase  modulation  (SPM)  effect  due  to  the  refractive  index  change  during  the  pulse 
modulation.  It  is  found  that  the  SPM  is  the  major  obstacle  that  limits  the  high  enough  power  mode-locked 
pulse  output  power.  Under  the  effect  of  SPM,  it  is  shown  that  the  high  power  operation  as  well  as  narrow 
pulse  generation  is  possible  by  increasing  the  absorber  saturation  effect. 
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4.0  Program  Conclusions  and  Suggestions  For  Future  Work 

The  work  carried  out  in  this  program  has  demonstrated  that  it  is  possible  to  fabricate  ultrashort  pulse,  wide 
spectral  bandwidth  semiconductor  modelocked  lasers.  We  also  analytically  showed  that  the  limiting 
property  of  these  lasers  that  prohibits  their  operation  with  many  (>  16)  phase  locked  longitudinal  modes  is 
the  occurrence  of  self  phase  modulation  (SPM)  in  the  dynamic  operation  of  the  laser.  The  analysis  also 
shows  that  for  a  laser  with  minimal  SPM,  i.e.  a  =  0,  it  will  be  possible  to  phase  lock  many  more  modes  and 
achieve  the  goal  set  out  in  the  conception  of  this  program.  This  zero-SPM  condition  in  principle  can  be 
achieved  on  a  round  trip  basis  by  careful  balancing  of  the  operating  points  of  the  various  sections  of  the 
CPM  laser.  It  is  unknown  at  present  whether  the  zero-SPM  condition  can  be  achieved  over  a  broad  spectral 
bandwidth  in  this  way.  More  extensive  modeling  will  be  necessary  to  determine  that  bandwidth  over  which 
this  condition  can  be  achieved. 

In  this  program  we  have  developed  the  epitaxial  tools  that  will  allow  us  to  tailor  the  energy  gap  in  various 
segments  of  the  laser  to  insure  that  the  round  trip  phase  modulation  is  minimized.  We  also  have  shown  that 
high  quality  devices  can  be  grown  using  these  techniques.  What  need  to  be  done  is  to  translate  accurate 
models  of  the  device  performance  to  real  physical  device  structures  that  would  allow  us  to  fabricate  these 
devices  and  test  the  ideas  experimentally.  Future  programs  should  undertake  this  activity. 

A  second  approach  that  may  bear  fruit  is  to  design  the  optical  feedback  structure  of  the  device  so  that  a  part 
of  the  SPM  is  cancelled  by  a  permanent  spectral  variation  in  the  phase  of  the  reflected  wave.  This  can  be 
accomplished  by  the  use  of  a  grating  or  photonic  crystal  mirror  in  which  the  reflectance  band  is  positioned 
so  that  the  high  dispersion  region  compensates  the  SPM  of  the  device. 

Finally,  it  may  be  possible  to  build  the  device  from  an  active  region  that  is  inherently  zero  SPM.  This  can  be 
accomplished  by  the  use  of  an  inhomogeneously  broadened  spectral  transition.  A  quantum  dot  active  region 
accomplishes  this  by  creating  an  ensemble  of  uncoupled  spectral  delta  function  emitters.  The  natural 
variation  of  the  dot  size  leads  to  this  spectral  dispersion  and  to  the  inhomogeneous  lineshape.  It  would  be 
interesting  to  analyze  the  consequences  of  the  use  of  QD  active  regions  in  this  context  and  to  explore  the 
fabrication  of  such  devices. 
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